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. Abstract. We show that the classical Bernstein polynomials Bj^i{f){x) on the interval [0, 1] 

P3 ' (and their higher dimensional generalizations on the simplex Em C M™) may be expressed 

^ . in terms of Bergman kernels for the Fubini-Study metric on CP™: BN{f){x) is obtained 

by applying the Toeplitz operator f{N^^Dg) to the Fubini-Study Bergman kernels. The 
expression generalizes immediately to any toric Kahler variety and Delzant polytope, and 
gives a novel definition of Bernstein 'polynomials' Bf^N (/) relative to any toric Kahler variety. 
, They uniformly approximate any continuous function / on the associated polytope P with 

all the properties of classical Bernstein polynomials. Upon integration over the polytope one 
obtains a complete asymptotic expansion for the Dedekind-Riemann sums J2aeNP fiw^ 
1^ ' of / e C°°(R™), of a type similar to the Euler-MacLaurin formulae. 



Introduction 

Our starting point is the observation that the classical Bernstein polynomials 



> 

QQ , aGN™:|a|<Af 

^ ■ on the m-simplex C M"* may be expressed in terms of the Bergman-Szego kernels 

g : Ilf,N^{z,w) for the Fubini-Study metric on CP"": Let e'^ denote the standard T™ = (S'^)"" 

C — ! action on C'" and and let Dg denote the linearization (or 'quantization') of its infinitesimal 

95 ; generators on //"(CP™, 0{N)). As will be shown in gH (see also g3]), 

where / G C^(M™). Here, Il^jiv denotes the Szego or Bergman kernel on powers 0{N) — >■ 

FS 

CP™ of the invariant hyperplane line bundle, f{N~^Dg) is defined by the spectral theorem 
and Hhps moment map corresponding to hps- Thus, the Bernstein polynomial B]\ff{x) 

is the Berezin lower symbol for the Toeplitz operator Uj^j f {N~^ D0)ll]\r , i.e. the value of its 
kernel on the diagonal. From this formula, many properties of Bernstein poljTiomials may 
be derived from properties of the Fubini-Study Bergman-Szego kernel. 

Furthermore, the formula ([2]) generalizes immediately to any polarized toric Kahler variety 
(L, M, u) and defines analogues B^n (/) (x) of Bernstein polynomials for any Delzant poytope 
P and any positively curved toric hermitian metric h on the invariant line bundle associated 
to P. We simply replace the Hermitian line bundle 0{1) —>■ CP"* with its Fubini-Study 
metric by any toric invariant Hermitian line bundle (L, h) — > {M,uj) (see Definition [2]). 
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The connection between Bernstein polynomials and Bergman-Szego kernels may be used 
to obtain asymptotic expansions of Bernstein polynomials as the degree N ^ oo; 

Theorem 0.1. Let [L, h) — > {M,uj) be a toric Hermitian invariant line bundle over a toric 
Kdhler variety with associated moment polytope P. Let f G C^{W^) and let BhN{f){x) 
denote its Bernstein polynomial approximation in the sense of Definition [H Then there 
exists a complete asymptotic expansion, 

BhHf){x) = fix) + Cifix)N-' + C2fix)N~^ + ■■■ + C„J{x)N-"' + 0(iV— i), 

in C°^{P), where Cj is a differential operator of order 2j depending only on curvature in- 
variants of the metric h; the expansion may be differentiated any number of times. 

In the case of classical Bernstein polynomials ([T]) (i.e. the interval or simplex) , this 
expansion has recently been derived by L. Hormander [Ho] by a different method (see [6]). The 
approach taken here is to use the Boutet de Monvel-Sjostrand approximations of Bergman- 
Szego kernels, with some simplifications in the case of toric hermitian metrics |BSj] ISTZ] . 
The operators Cj are computable from the coefficients of the asymptotic expansion of the 
Bergman-Szego kernel IlfiN{z,z) on the diagonal in |Z2^ ILuj . It should be noted that for 
general toric Hermitian line bundles, the Bernstein 'polynomials' are not quite polynomials 
in the usual sense, although they are algebro-geometric objects in the sense of [Dl|T]; see §2] 
for further discussion. 

As defined in ([2]) and in Definition [21 the Bernstein polynomials are quotients 

Bh^{f){x) = . . -li ^^ (3) 

of a numerator polynomial f/f^N f (x) by the denominator IIi^n{z,z) with ^h{z) = x. Here, /x^ 
is the moment map associated to the Kahler form Uh associated to h. The numerator poly- 
nomials also admit complete asymptotic expansions, and indeed the Bernstein polynomial 
expansions are derived from the numerator expansion and from the asymptotic expansion of 
the denominator. Hence, Theorem lU.ll follows from: 

Theorem 0.2. With the same assumptions as above, there exist differential operators Afj, 
such that 

K-{fKx) ~ N^fix) + N^-'j\r,f{x) + ■ ■ ■ , 

where the operators Mj are computable from the Bergman kernel expansion for H/jJv(2;, z). 

Theorem 10.21 has an application to Dedekind-Riemann sums over lattice points in dilates 
of the polytope P, i.e. sums of the form 

Upon integration of M^n f{x) over P one obtains: 

Corollary 1. Let f G C^(M™). Then there exist differential operators Sj, such that 
Yl /(^) ~ / f(^)dx + ^— / f{x)da + iV"-2 / S2f{x)dx + ■ ■ ■ , 

aeNP ^ JP ^ JdP JP 

where a is the Leray measure on dP corresponding to the affine defining functions ir{x) = 
{x, Vr) of the boundary facts (cf. \2B^) . That is, on the rth facet of dP, dir A da = dx. 
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Exact and asymptotic formulae for J2aeNP f(jf) have been previously proved for special 
/ using the generalized Euler-MacLaurin formulae of Khovanskii-Pukhlikov, Brion-Vergne, 
Guillemin-Sternberg and others (cf. [Gl IGSl IGSWt IKSWj ). For purposes of comparison, 
Theorem 4.2 of ^ states that for / e Co°°(K"), 

E /(^) - 1 E E -^(^^m) / f^^^A U.0 + O(iV-), (4) 

where the sums involve various data associated to the polytope P and where Ph is a parallel 
dilate of P. We refer to [GSJ for the notation. The two term expansion given in Corollary 
[1] was stated in [Sz] . It is straightforward to generalize the formula and proof to the case 
where / is a symbol as in [GSWj . and to obtain remainder estimates in the expansion. 

A significant difference between the Euler-MacLaurin and the Bernstein methods for ob- 
taining expansions of Dedekind-Riemann sums J2aeN'^-\a\<N fiji) ^^^^ Bernstein ap- 
proaches uses an arbitrary toric Kahler metric while the Euler-MacLaurin approach is metric 
independent. This reflects the fact that the Bernstein approach is to integrate the pointwise 
expansion of Theorem 10. 2^ which depends on the metric h. The metric independence of 
the expansion in Corollary [1] is equivalent to a sequence of integration by parts identities 
involving curvature invariants. For instance, we obtain the second term in the expansion in 
^ by using an integration by parts identity on polytopes due to Donaldson |D2] ; see also 
^ for the simplest case. Conversely, comparision of the metric expansion in Theorem 10.21 
and the Euler-MacLaurin expansion in (jl]) gives another proof of this identity, and generates 
further identities in the lower order terms for any choice of toric hermitian metric. 

The connection between Bernstein polynomials, Bergman kernels and Berezin symbols 
appears to be new, and one of the principal motivations of this article is simply to point out 
the toric geometry underlying the classical Bernstein polynomials. But a further motivation 
is that the generalized Bernstein polynomials should be useful in the program of Yau-Tian- 
Donaldson of making algebro-geometric (i.e. polynomial) approximations to transcendental 
geometric objects on Kahler varieties (cf. [Pit [T]). For instance, in [SoZt ISoZ2] what we 
recognize in this article as Bernstein polynomials were used to approximate geodesic rays in 
(see also |PSj ). However, the function / in that paper also depended on in a subtle 
way and so the polynomials were much more complicated than the Bernstein polynomials 
of this article. The article [Ho] also concerns relations between Bernstein polynomials and 
Bergman kernels, but mainly for the opposite purpose of deriving Bergman kernel expansions 
on Reinhardt domains from classical Bernstein polynomial expansions on the simplex. The 
exposition in ^ was influenced by its analysis of Bernstein polynomials. It also draws on 
some of the analysis of [SoZj . 

In addition to the Bergman-toric generalization of Bernstein polynomials, there also exists 
a probabilistic generalization of Bernstein polynomial which replaces (^) by the weighted 
number of lattice paths from to a with steps in the polytope P. This definition also 
coincides with the canonical one in the case of the Fubini-Study metrics on CP™ but in 
general gives a different class of polynomials defined on the simplex of probability measures 
on {!,..., m}. In the case of the simplex = P, both spaces are the same, but in general 
they are not. The relevant analysis could be obtained form [TZ]; we will not discuss these 
generalizations here. 



4 



STEVE ZELDITCH 



We would like to thank H. Hezari for a careful reading of the article and for pointing out 
some notational inconsistencies and misprints in an earlier version. 

1. Fubini-Study and classical Bernstein polynomials 

Let us begin by explaining in more detail the Bernstein-Bergman connection for the Fubini- 
Study metric in one complex dimension. We recall that Bernstein polynomials of one variable 
give canonical uniform polynomial approximations to continuous functions / G C([0, 1]): 

BN{f)ix) = J2( ■ /(^y(l-x)^-^. (5) 

They have the special feature that they simultaneously uniformly approximate all derivatives 
of / if / e C^ i.e. f^\x) (cf. [Lj), and if / G there exists a complete 

asymptotic expansion ( |Hoj ) 

oo , 

^J.=o 

for certain polynomial differential operators L^[x, 
In this case, 

BnU) = Wn^Nif) (cf. Theorem [021), and also 



jJJo ' ' \jj {N + l)\ N+1 

Hence, (jSD implies that 



f^'AfM){x)dx = E;Lo/(i) 

= (AT + 1) (/^^ f{x)dx + ^ J^\x - x')f"{x)dx + ■ ■ 

= {N+1) f{x)dx + ^ (/(I) - /(O) - 2 J^' f{x)dx) + 

= Nj^f{x)dx + \uiX)-fm+o{j^). 

(7) 

We included the routine details to point out that obtaining the first two terms of the Euler- 
MacLaurin Riemann sum expansion in Theorem [1] required two integrations by parts and 
cancellations of f{x)dx in the constant term between the subleading term of the dimension 

(Riemann-Roch) polynomial (A^-|-l) term and in the Lif(x)dx term. Similar cancellations 
occur in the general case (see the proof of Theorem [T|). 

We now relate the Bernstein polynomials B^lf) on [0, 1] to the Bergman kernel for the 
Fubini-Study metric on CP^. The discussion is almost the same for the m-simplex C M™' 
and the Bergman kernel for the Fubini metric on CP™, so we carry it out in all dimensions. 
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We first need to recall some standard facts about the Bergman or Szego kernels for the 
Fubini-Study metric. 

By the m-simplex we mean the convex set = {(a^i, • • • ,Xm) ^ '■ \\x\ \ := < 
1}. We denote its dilate by G M by NJ^m- As discussed in |STZ] and elsewhere (see 
[STZj for references), the space Voly{NT^rn) of polynomials with exponents a G NT^m can 
be identified with the space of degree- iV homogeneous holomorphic polynomials in m + 1 
variables by identifying the (non-homogeneous) polynomial 

/(^i, . . . , z^) = C"^" = ■■■zT) 

\a\<N 

with the homogeneous polynomial 

-^(Co, • • • 5 Cm) = ^ CaC,Q '"'Cr^ ■ ■ ■ Cm™ " 
|Q:|<Af 

The space Voly{NT,m) has a natural inner product, 

{f,9) = ^[ FGdu, (8) 

ml Js2m + 1 

This inner product is equivalent to viewing f,g as a holomorphic sections of the iVth power 
0{N) of the hyperplane line bundle C(l) — > CP'" dual to the tautological line bundle. The 
line bundle 0{1) carries a natural metric hps given by 

I|s|Ups(M) = , ^ = K, ...,wje c™+^ , (9) 

for s G C^+i* = if0(CP'",O(l)), where |wp = EJLo l^jf and [w] G CP'" denotes the 
complex line through w. The Kahler form on CP"* is the Fubini-Study form 

O^FS = ^e,,, = y^dd\0g\w\\ (10) 

The natural Fubini-Study inner product on sections is then 

In an affine chart and local frame e, sections have the form fe where / is a polynomial and 
the inner product takes the explicit form 

{fr9) = -J f,9eroiy{Nj:^). (ii) 

mi jQm [i -\- \\z\\ ) 

Both versions of the inner product generalize to any holomorphic line bundle. 

A basis for Voly{NT,m) is given by the monomials Xa{z) = z^^ ■ ■ ■ z"^ , \a\ < N. The 
monomials {Xa} are orthogonal but not normalized. Their norms given by the inner 
product ([8]) are: 

1 

(12) 



\Xa\\ 



{N-\a\)\a,\---aJ. 
(N + my. 
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Thus, an orthonormal basis for Voly{N'Em) is given by the monomials 

(iV + m)! 



\\Xa\ 

where 



1 



{N-\a\)\ai\---aJ. 



. (N + m)\ fN\ , , 



(14) 



a J {N — \a\)\ai\ ■ ■ ■ am^. 
We let Xa '■ -S"^"*"*"^ — > C denote the homogenization of Xa- 

x^(x)=xriv (15) 

The Bergman or Szego kernel II^jjv for the Fubini-Study metric is the orthogonal projection 

FS 

to the space H'^{CF"^,0{N)) of holomorphic sections with respect to the inner produced 
induced by hps, which hfts to the orthogonal projection to Voly{NT). It is thus given by 

Il^{x,y) = J2 ^Xa(^)Xc:M= ^^^^^(^,^)"^ , (16) 

|a|<Ar 

for x,y & S^'^'^^. In particular, on the diagonal we have {x,x) = 1 and 

n,(x,.) = <^. (17) 

In terms of the standard local affine frame on C"*, we have Xa (^) = ^i+\\l"\2y/2 : and hence 

^ , , (iV + m)! EH<iv(^)^"^" 

Wc now have the ingredients to identify Bernstein polynomials for the simplex NT^m in 
terms of the Fubini-Study Bergman-Szego kernel. The Kahler potential of the Fubini-Study 
metric is (pps — log(l -|- \ \z\\^) where = J2j kjl^' i^^ moment map is 

U 12 

/^ftFs(^) = ( 

The Fubini-Study symplectic potential is the convex function on given by the Legendre 
transform of (pps in logarithm coordinates, 

m 

^oi^) = X^^jlog^j + (1 - lkll)log(l - \\x\\) 
i=i 

where | |a;| | = Xj. A simple calculation shows that the Bernstein terms may be expressed 
in terms of the symplectic potential as 

U"(l - ||x||)^-H - „ - ^ . (19) 

It follows that 







2 


Z-m 


2 


1 + 1 


z\ 


r---'(i+i 


z\ 


2 



BN(f)(x) = n-rcZiT Elo m) ^ - - ^hU^)- (20) 



"■FS 
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On the other hand, one can also express the Bergman- Szego kernel in terms of the symplectic 
potential at the points {e^^z, z) as 

n (^^e^ A - iv(„o(.)+(t-.,v.oW>) 



(21) 

^h-Jz,z)J:Lo Oe^^"a;"(l-||a:||)^-H. 



Indeed, comparing (fTSi) and (12T|) . we see that the two expressions for the Bergman-Szego 
kernel agree as long as 

|^a|2g-7Vlog{l+||.|n ^ ^;v(no(x)+(t-x,V.o(.)»^ ^hcn fX,^^{z) = X, (22) 

and this follows from the pair of identities, 

|^a|2 _ ^{a,Vuo{x))^ j^g^^ ^ |^|2^) _ VUo(a;)) - %(x) whcn /i/,^5(2) = X. 

On the open orbit, we may use logarithmic coordinates z = ePl"^^"^^ . Then p = '\7uo{x) and 
the identities are equivalent to the fact that the Kahler potential and symplectic potential 
are Legendre transforms of each other. Since both sides of (PTi) are continuous, the equality 
extends to all of M and P. 



Applying the operator f{jf-) just replaces e*^° by f{-^)- Then, dividing by Uj^ 



N (z,z) 

FS 



gives ( !20|) and (12]) . Together with the formulae above for norms of monomials and the Szego 
kernel in dimension m, the formula (tH) also reduces to ( |20D . 



2. Definition of the generalized Bernstein polynomials 

We now generalize the definition of Bernstein polynomial to any polarized toric Kahler 
variety, and generalize the calculations of the previous section. 

We recall that a toric Kahler manifold is a Kahler manifold (M, J, u) on which the complex 
torus (C*)'" acts holomorphically with an open orbit M". We assume that M is projective 
and that P is a Delzant polytope, i.e. a convex integral polytope in with the property 
that each vertex is contained in exactly m facets, and the normals to the m facets at each 
vertex form a Z-basis for a lattice F C so that T™' = M^/F is the torus acting on Mp. 
The convex polytope P is defined by a set of inequalities of 

{x,Vr) > Xr, r = l,...,d, (23) 

where Vr is a primitive element of the lattice and inward-pointing normal to the r-th (n — l)- 
dimensional face of P. 

We denote by T'" = (S^)"^ the real torus underlying (C*)'". By a toric Kahler metric we 
mean a Kahler metric uj invariant under T*". We assume that is a de Rham representative 
of the Chern class ci(L) e if^(M, M) of an invariant holomorphic line bundle L — ^ M. We 
let h denote the Hermitian metric on L inducing the Chern connection with curvature (1, 1) 
form uJh = oj. Here, given a Hermitian metric 

^h = -^dd\og\\eL\\l, (24) 

where cl denotes a local holomorphic frame (i.e. a nonvanishing section) of L over an open 
set U C M, and WeiWh = h{eL, eiY^"^ denotes the /i-norm of e^- We often write cu for tOh 
when the metric is fixed. 
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Now fix a basepoint mo on the open orbit and identify M° = (C*)"^, endowing M° witli 
tlie logaritlimic coordinates 

z = eP/^+'^ G (C*)'", p, ^ G 

Over tlie open orbit, u lias a Kahler potential, i.e. u = —2idd(f{z). The associated Hermitian 
metric then has the form h = e"*^. Invariance under the real torus action implies that (p 
only depends on the p-variables, hence, 

uj = - y dzj A dzk- 

2 ^ dpkPj 

We sometimes subscript u to indicate the associated hermitian metric or Kahler potential, 
e.g. uj = ujh = uj^. By a slight abuse of notation, we denote the Kahler potential in the 
logarithmic coordinates by ip{p). Positivity of u implies that (f is strictly convex of p G M". 

The real torus T*" acts on (M, a;) in a Hamiltonian fashion with respect to u, and its 
moment map /x<^ = with respect to = ojh is defined by 

= Vp¥^(pi,...,p™), {z = eP'^+''). (25) 

The symplectic potential u^p associated to the Kahler potential is defined to be the Legendre- 
dual of if, defined as follows: for a; G P there is a unique p such that /i^(e''/2) = V p^p = x. 
Then the Legendre transform is defined to be the convex function 

u^{x) = (x, p) - ^(p), ^ (26) 

on P. 

There exists a 'canonical' Kahler metric and symplectic potential, defined as follows: Let 
Ir : R" — > R be the affine functions, 

ir{x) = {x, Vr) — Xr- 

Then the canonical symplectic potential is defined by 

= ^4(x)log4(x), (27) 

k 

which in turn corresponds to a canoncial Kahler potential [Gl|X]. Every symplectic potential 
has the same singularities on the boundary dP as the canonical symplectic potential. 

We denote by = V^u^ the Hessian of the symplectic potential. It has simple poles on 
dP. We also denote by H^^p) = VpV9(e''/2) the Hessian of the Kahler potential on the open 
orbit in p coordinates. By Legendre duality, 

Hpip)=G-\x), fxpie^/')=x. (28) 

We now let (L, h) ^ M denote the invariant Hermitian line bundle with curvature Uh = uj- 
A natural basis of the space of holomorphic sections H^{M, L^) associated to the A^th power 
of L ^ M corresponds to monomials 2;° where a is a lattice point in the A^th dilate of the 
polytope, a G NP fl Z"^. The hermitian metric h on L induces inner products Hilb]\f{h) on 
H%M,L^), defined by 

(Sl,S2)/jiV= / {Si{z),S2{z))hN-^. 
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The monomials are orthogonal with respect to any such toric inner product and have the 
norm-squares 

Qh^{a) = [ |z"|2e-^'^(^)rfV;(^), (29) 
where dV^ = {iddip)'^ /m\. In terms of the symplectic potential, 

Qf^N{a) = J e^("^(")+<t-.v-^(-))rfa;. (30) 

The Bergman-Szego kernels for this hermitian metric are the orthogonal projections with 
respect to HilhN{h) to H^{M, L^). If we denote the sections corresponding to the monomials 
by Sa then, 

UhN{z,w)= > — — . 

The following definition generalizes the formula of ([I]) to any toric Kahler manifold. 

Definition: Let / G C{P). The A^th normalized (Bergman-)Bernstein polynomial approxi- 
mation to / with respect to the hermitian metric /i on L — ^ M is defined by 

BhNf{x) = u^^^fh^fix), where 

As in the classical case, Bernstein polynomials are closely related to certain probability 
measures on P. We define 

;.^:=V^4^5., (31) 
.t^pn,(z,z) 

where VhN{a,z) denote the Fourier coefficients of the Bergman kernel with respect to the 

Vf,N{a,z):= ' ' . ■ (32) 
Proposition 2.1. Let f e C{P) and let x = iJ„p{z) and let h = e~'^. Then, 

aeNP J I AT/ nNiz,z) ' 

Proof. The first two equalities are obvious from the definition. The third equality generalizes 
the identity ([22]): 

|^a|2g-7V^(.) ^ g7V(«^(x)+(f-x,log^-i{x)»^ ^j^g^ ^ ^_ ^33^ 

As in the case of the Fubini-Study metric, the identity splits into two identities on the open 
orbit, 

|^a|2 ^ g(a,p)^ ^-N^{z) ^ ^N{u^(x)-{x, log fi-^(x))) ^ ^3^^ 
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The first follows from the fact that 

Vx-M^(x) = log^-^ix) = p, (35) 

since by ( 126|) . VxU^{x) = p + {x, V xP) — {^f{p), V^p) = p, as Vip{p) = x. The second then 
follows from the fact that <f{p) and 'U^(x) are Legendre duals. The identity of the Proposition 
then extends by continuity to the closure. 

□ 

As a simple corollary, we obtain one of the standard properties of Bernstein polynomials. 

Corollary 2. Let f e C{P). Then minp/ < BN{f){x) < maxp/. 

Let us calculate explicitly the numerator polynomials for the canonical symplectic potential 
(127|) or Kahler form. We have, 

k k 

Hence, 

^N~^' '^«o(a;)) = 2_^{—-x,Vk)\ogik + {j^-x,v), 

k 

and 

where in the last line we use that ik{x) — {x,Vk) = —Xk- Hence, the numerator of the 
canonical Bernstein polynomial may be rewritten as 

K-fix) = EaeNP /(# )Q^e<"-^-'^)H,(4(x))-^^^+<'^'^^>, (36) 
which closely resembles the classical cases (where also v = 0). Here, 

Q^Nja) = J e^("oW+(f-x,Vno(x)>^3, 

is the norming constant with respect to the canonical symplectic potential. 
In general, the symplectic potential has the form 

u^{x) = uq{x) + g^{x) = ^4(a;) log4(x) + fi'^(x), (37) 

k 

where g^^ G C°°{P) is smooth up the boundary [G ] R l [D2] . Hence the a term gets multiplied 
by the additional factor 

gAf(g^(x)+(f-x,Vg^(x)>)_ 

In Definition 10.11 Bernstein polynomials were normalized by dividing by H^jv(2;, z). It 
follows that i?/jjv(l) = 1 as for classical Bernstein polynomials. However, in the classical 
cases, IlhN{z,z) is constant so the normalization of the polynomial is still a polynomial. In 
general, however, U^n [z, z) is not constant on the diagonal and therefore the quotient is rarely 
a polynomial in the usual sense. In special cases, Bernstein polynomials are polynomials in 
the variables x = p{z), but this depends on the properties of the moment map of the Kahler 
metric. 
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One might prefer to normalize by dividing the numerator polynomial by the dimension 
polynomial d]sr = dim H^{M, L^), 

B,Nf{x) = i^KNf{x). (38) 

For the canonical metric, one would then have the canonical Bernstein polynomials 

Bnr^m = iE„.;vp/(f )o;i(^e<--^^'')n,(4(x))-^^^+<"'-), (39) 

which are visibly polynomials when v = 0. However, this is essentially an aesthetic decision 
based on how seriously one wants to take the term 'polynomial'. Either definition has the 
same value in terms of making approximations. Our view is that the term 'polynomial' should 
have the same general sense in the Kahler context as 'algebro-geo metric' approximations do 
in the Yau-Tian-Donaldson program. 

To our knowledge the only previously studied cases are the Bernstein polynomials for the 
simplex ([1]) or cube, 

BnUHx) = E /(|, . . . , ^)nti (f - x,)''-'^. 

0<h,...,ia<N ^ 

Here, {xi, . . . ,Xd) G [0,1]'^. The classical Bernstein polynomials B]\f{f) are distinguished 
among other polynomial approximations by simultaneously approximating the derivatives 
and also by preserving certain shape and convexity properties (at least, in dimension one). 
It might be interesting to explore the shape preserving properties of general Bernstein poly- 
nomials. We also note that Bernstein polynomials admit holomorphic extensions to complex 
neighborhoods of the polytope P when the symplectic potential function is real analytic. 

3. Bernstein polynomials, Toeplitz operators and Berezin symbols 

In this section, we prove formula ([2]) and also establish some basic properties of Bernstein 
polynomials. 

The proof of is simply a matter of unwinding the definitions. The Bergman kernel is 
a section of the bundle {L^) (L^)* ^ M x M. It is simpler to deal with scalar kernels, 
and so we lift the Bergman kernel to a kernel IIiv(x,y) on the unit circle bundle X — > M 
with respect to h in the dual line bundle L*. In other words, X = dDf^ is the boundary 
of the unit disc bundle with respect to h in the dual line bundle L*. We use local product 
coordinates x = {z,9) G M x on X where x = e^^ \\}(z)\\ terms of a local holomorphic 
frame e{z) for L. When working on M we tacitly use the representative of H^n relative to 
the frame e{z)^ of . For the sake of brevity, we will not review the definitions but refer 
to |STZj for the relevant background. 

The space H^{M, L^) is naturally isomorphic to the space Hff{X) of CR holomorphic 
functions transforming by e*^^ under the action of the circle bundle X ^ M. We denote 
by s — > s the lift of a section to an equivariant CR function and by II^jv(x, y) the lifted 
Szego kernel, i.e. the orthogonal projection from L'^{X) — >■ if^(X). The monomial sections 
Sa which equal z°' on the open orbit lift to equivariant functions Sa on X. 

By the standard linearization of geometric quantization (reviewed in this context in |STZj ). 
the T*" action lifts to X as contact transformations of the Chern connection form associated 
to h. For the sake of completeness, let us recall the lift of the torus action to Hjf{X), and its 
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linearization on H^{M, L^): The generators ^ of the T"* action on M hft to contact vector 
fields Si, . . . on X. There is a natural contact 1-form a on X defined by the Hermitian 
connection 1-form, which satisfies da = 7c*uj. The horizontal lifts of the Hamilton vector 
fields are then defined by 

7r4 = e„ «(e') = o, 

and the contact vector fields Ej are given by: 

I, =i^ + 2vrz(/i o vr,^;)- = + 2vrz(/i o vr),- -, 

where is the moment map corresponding to h, and where ^* G is the element of the 
Lie algebra of T™" which acts as on M. 

It follows that the vector fields act as differential operators on the CR Hardy spaces, 
: H%{X) ^ H%{X) satisfying 

(S,5)(C) = ^^(e^^-C)U=o, SeC°°{X^p). (40) 

Furthermore, the generator of the 5*^ action acts on these spaces and 

^■.nUX'p)^'Hl{X'p) , ll-s^ = Nsj, for s^enliXf,) . (41) 



Since by (1401) . the operators act by translating functions by the T"^ action lifted to X, 
we henceforth denote iS^ by Dq^. Then for 1 < j < m, the lifted monomials Xa ^ H'^{X) 
are joint eigenfunctions of these commuting operators, 

De.Xa = ajXa, Va E NP. 

The dilation P NP is best viewed in terms of constructing a conic set of eigenvalues in 
one higher dimension by adding the operator 

W = f^-E^^r (42) 

i=i 

The monomials xs are then the joint eigenfunctions of these (m + 1) commuting operators 
and we define the 'homogenization' NP C Z™'"'"^ of the lattice points in the polytope NP to 
be the set of all lattice point of the form 

= a:= {ai,...,am,N -\a\), a = {ai, . . . ,am) e NP nZ"^, (43) 

Given / G C°°(M'"), we now define f{De) on L^{X) by the spectral theorem for m com- 
muting operators, i.e. 

f{De)=[ /(Oe^<«'^*>rfe, where {^,De) = J2^^De^- 



We then have 



fiN-'De)s^ = f{^)s^ (44) 



Since UhN{z,w) = ZlagArp ■5Q(^)'5a(t&), we have 

V fi 

'N 



a&NP 

f{N''De)nnN{e''z,w)= ^ f[^)s^[z)J:^). (45) 



aeNP 
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It follows that 

fiN-'De)Il,Mie^'z,w)U=^ = /(^) (46) 

aeNP 

The right hand side is constant along the orbits of the action and may be identified with 
a function of z G M. On M we have |so(5)|^ = ||sa(2;)||^jv and by Proposition 12. II we obtain 
the definition of the numerator polynomials when we substitute z = yU^^(x). Equivalently, 

Ar,.(/)(x) = (fl,NfiN-'De)h,M) (e^^^,z)|,^o^,^^-i(,), (47) 

where the right side is the Berezin symbol of the Toeplitz operator Ylf^N f{N~^Dg)IlhN . We 
then divide by n^jv(2, z) to obtain the Bernstein polynomials. 

4. Proof of Theorems 10.11 and 10.21 

We now use the Boutet de Monvel - Sjostrand parametrix |BSj[ |BerSj[ |BBSj| to obtain a 
complete asymptotic expansion for the Bernstein polynomials from (H71) . There now exist 
many expositions of the construction and properties of this parametrix, so we will only briefly 
recall the essential elements in the case of toric varieties [SoZUSTZ] . We also use the notation 
X, y for points of X, hoping that no confusion with coordinates on P will occur. 

We first recall that, on the diagonal, the Bergman-Szego kernel has a complete asymptotic 
expansion, 

UMz, z) = J2\\S^{z)\\l^ = ^ + «i W^'' + «2(^)iV-2 + ■■■], (48) 

i=0 

for certain smooth coefficients aj{z). In fact, 

{a ~ -S 
al = IaS + ^(|i?|2 - 4|fficp + 3S^) ^^^^ 

where R, Ric and S denotes the curvature tensor, the Ricci curvature and the scalar curvature 
of ujh, respectively, and A denotes the Laplace operator of (M, cu/j).; see [Z2l|Lul BSj BBSj] . 
Off the diagonal we have the following expansion: 

Proposition 4.1. For any C°° positive hermitian line bundle {L,h), there exists a semi- 
classical amplitude in the parameter , siy{z,w) N"^so{z,w) + N"^~^si{z,w) + - ■ ■ , such 
that 

U^,n{z, w) = e^(^(^'-)-i(^(^)+^('")))s^(z, w) + 0{N-°^), 

where (p is a smooth local Kdhler potential for h, and where ip{z,w) is the almost- analytic 
extension of (p{z) = ip{z,z). 

Since the local Kahler potentials (e.g. the Kahler potential on the open orbit) are invariant 
under the T*" action, they can be expressed in the form F(|2;p) where F e C°°(M). We 
denote by F{z ■ w) the almost analytic extension of F. Thus, we have: 

Proposition 4.2. For any hermitian toric positive line bundle over a toric variety, the 
Szego kernel for the metrics have the asymptotic expansions in a local frame on M , 

Ii^N{z,w) ~ e^(^(^-'^)-^(^(ll^ll')+^(ll'"ll'»)A;v(^,^) modN-^, 
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where An{z, w) ~ N'"^ ^1 + "■'^^^'^^ + ■ ■ ■ j is a semi-classical symbol of order m. 

We now prove Theorems 10.11 and I0.2[ 
Proof. We then apply the geometric quantizations of the torus action to get, by Definition 

m 

7^ — r~{ • 

aeNPnZ"^ ^ ' 

By P7I) . we obtain Af^Nf^x) by integrating the right side against /(O- We note that in 

general e*^^'^ ^^®V(c*^u^)|e=o = ^(e*^^"*"^^w)|e=o = 'ip{e'^Nw) Performing the same transfor- 
mation on the parametrix gives, 

ATMfKx)- [ /(Oe^(^(^'""^'^'^^-^(l^l^)M^(.,e4.)d^, (50) 



where ~ means that the difference is a function which decays rapidly in along with its 
derivatives. Such a remainder may be neglected if we only consider expansions modulo 
rapidly decaying functions of A^. 
We have, 

= iN-' /o (V^F(e^*^"'«+^), e)dt 



(51) 



where 



R^{tN,z) := /'(t- l)2vJ(F(e'*f+0)(^OV3!, (52) 
Jo 

and where = V^F(|zp) = V^(/?(e^) is the Hessian in the notation (1251) . Hence, fl5U]) takes 
the form 

BN{f){x) ~ 4„/(Oe'^^(^)'^>e(^^)"'<^^€'«>+^"'«3(«,^,.)^^(^^e4z,0,Ar)c/^ (53) 

and by Taylor expanding the factor e'-*^^ ^{HzS.,0+^ ^R-i{d,N,z) obtains an amplitude A^^ 
such that 

MhN{f){x) ~ ^,„/(Oe^<'^(^)'«>i7v(^,e4^,O,iV)d0. (54) 
The amplitude Aj^ has an expansion of the form. 
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for various smooth coefficients aj{z); the ffist one is constant. If we divide by U^nIz, z) we 
cancel the constant and by expanding the denominator we obtain, 

KmU){x) ~ iV™/(/i(^)) + iV--i {i^HH^D^, D,)f{fi{z)) + a^{z, z)f{fi{z))) + 0(iV™-2), 

(55) 

Since fj,{z) = x we obtain Theorem 10. 2[ Dividing by IlhN{z,z) and using f HHl) completes the 
proof of Theorem 10.11 □ 

It is difficult (but possible) to calculate the coefficients in explicit geometric terms by this 
method. In the next section, we will reduce the calculation to the known calculation of 
Bergman kernel expansion coefficients. 

4.1. Proof of Corollary [H To prove the Corollary, we integrate the expansion (155!) over 
P to obtain 

fpAfhN{f){x)dx = jpf{x)dx + 0{N^-^), (56) 

where the lower order terms could be computed from the expansion. But we postpone their 
evaluation until the next section. 

5. Bergman kernel expansion and geometric expressions for the Bernstein 

expansion of theorem [2] 

In this section, we give a second proof of the convergence of Bernstein polynomials which 
is based on one of their essential features: the localization of the sum over G P fl 
around the image of z under the moment map. This is well-known and various expositions 
can be found in |Ho[ [K | [L] : see also [D] Lemma 6.3.5). This approach reduces the calcinate 
the lower order terms in the Bernstein polynomial expansion in terms of the Bergman kernel 
expansion in [?, iLuj and elsewhere. 



The relevant Localization Lemma was proved in [SoZj. We use a notation similar to [Ho]. 
Lemma 5.1. (Localization of Sums) |SoZ] Let f G C{P). Then, there exists C > so that 

Hence, it is natural to Taylor expand / around /i/i(;z) to obtain 

/(a/iV)= /('^H/^.(eO)(^-MeO)7^! + ^M(/,e^^), 

v<2M 

where Rm is the Mth order Taylor remainder. We then have, 

( 



2 

N_ 

(a) 



(57) 



+7^(M, N,z), 

where the remainder is obtained by summing Ruifi e'', ^) in the variable ^. 
To prove the main result, we need to study the special functions 
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Proposition 5.2. Uniformly for z e M we have: 

I-^{z) = 0{N'^-''/\logNy). (59) 
Proof. The Localization lemma implies that 



ae7VPnZ'":|f-Mh{^)l<^^ 

In the domain of summation we then have, 



and this implies the statement. 



□ 



We can explicitly evaluate these functions by relating them to derivatives of the Bergman- 
Szego kernels. The following Lemma was also used in [SoZ] . We employ a tensor product 
notation - ^^[eP'^))f for - ^^h{ePl^)i)Q - ^^h{ep'^)j). In the following, we implicitly 
assume that z lies in the open orbit and express it as 2; = e^/^"*"*^. Similar formula hold at the 
boundary as well where the vector fields are replaced by derivatives in affine coordinates. 

For the sake of brevity we refer to |SoZ] for the modifications to the formulae around the 
boundary. 

Proposition 5.3. We have: 

Proof. To prove (1), we differentiate (l25l) to obtain 

To prove (2), we take a second derivative of (1) in p to get 



We now evaluate these functions geometrically: 

Proposition 5.4. We have: 

(1) ij^^iz) = C^N^-^VS{z) + 0(iV"-3); 

(2) = AT'^^iy^v? + N'^-^S{z)Vlip + 0(Ar"-3). 



□ 
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Proof. From fHSl) it follows that 

VpfiH{z)Ilj,M{z, z) = N^Vfih + N^'^CmS{z)V^ih + 0(iV— 2); 
We also use that Vfih{eP/^)) = VV- 



□ 



To complete the second proof of Theorem 10.11 it suffices to observe that the remainder in 
flFr|) after expanding to order M is 0{N''^~^^^'^{logN)^'^), which follows from the fact that 
n{M,N,z) < CfN'^rJ^z). Therefore 

X,.{f)if^{z)) = fifi{z))n,.iz,z) 

= N^f{ii{z)) + iV— 1 {f{ii{z))S{z) + Viih ■ VV(/u(^))-) + 0(iV--3/2(logiV)3). 

(60) 

6. DEDEKIND-RIEMANN sums over LATTICE POINTS: PROOF OF COROLLARY [1] 

As noted in §4.11 the existence of an asymptotic expansion for the Riemann sums follows 
immediately from theorem l0.2[ However, it is an expansion in terms of integrals of curvature 
invariants against derivatives of / over P. The purpose of this section is to prove that the 
first two terms can be put in the form stated in Corollary [T], and thus to clarify the relation 
between the Bernstein and Euler-MacLaurin approaches to lattice point sums. 

We begin the calculation by integrating (I6T!) over M with respect to ^ and recalling that 
the pushforward to P of this volume form under the moment map /i^ is Lebesgue measure 
dx on P. Also, lifiN^z^z) is constant on T™- orbits, so n/jjv(/i~^(x), is well-defined 
although the inverse image is an orbit. The same is true for geometric functions such as the 
scalar curvature. We also recall that = dim if°(M, L^)). Then by Proposition 15.41 only 
the zeroth and second order terms of the Taylor expansion of / contribute to the A^~^ term 
of the Riemann sum expansion, and we have 



+ \ Y.\p\=2 Ip D^f{x)liAx)dx + 0(Ar-3/2(iog Nf) 
= ^ jpf{x)dx+^ j^\f{x)S{^~\x)) 

+\{Vp^lH{^i-\x),Vlf{x))dx + 0{N^-^'\\ogNf). 
Here, (V/i/i, V2/(/i(z))) denotes the Hilbert-Schmidt inner product of the tensors. 



(61) 
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By Legendre duality, the Hessians of the Kahler potential and symplectic potentials are 
inverses, i.e. 

Vp^lH{^l'\x) = {v^u^{x))-\ (62) 

Hence, 

{VpiiH{li-\x),Vlf{x))dx= [ y^u^'f.kdx. (63) 

Further, we recall (cf. |D2t R] ) that the scalar curvature of a toric Kahler metric is given in 
terms of the symplectic potential by 

j,k 

where u^^, ^ ^ j, k < n are the entries of the inverse of the matrix V^u^. See [m] (3.1.4). 
We now use the following integration by parts formula due to Donaldson: 

Lemma 6.1. (^ |D2j . Lemma 3.3.5) For any symplectic potential u^^ and f G , ^^j^u^ fjk ^ 
L^{P) and 

I Y.<f^^^= I Y.K')^3kfdx+ I fda, 

JP JP JdP 

where da is the measure defined in CorollaryUl 
Combining Lemma [6. II and flMl) we obtain 

^-f{x)S{fi-\x)) + ^(V,/i,(/i-^(x), Vlf{x))dx =\j^^ fda, 

proving that the two term expansion in Corollary [T] is correct. 
Remarks: 

(i) We note that in |D2j Lemma 3.3.5, the boundary term is given the — sign. However, the 
measure da was only defined there (page 307) up to sign. The sign of this term is universal 
and by comparing with the one-dimensional case, we see that it is positive. 

(ii) To connect this calculation to the classical one- dimensional case ([7]), and perhaps clarify 
the notation, we note that its N"^~^ (with m = 1), 

f(x)dx+- / (x — x'^)f"{x)dx, 
2 Jo 

may be expressed in terms of the Fubini-Study Kahler potential and moment map as 

- -^fiFs{lJ'~^{x))f"{x)dx, x = n{eP^^), 

"P 

since 

^Me"/') = log(l + eO, ^VFs{eP'^) = f^rsie"/') = = x, 

dp I + eP 

and 

d"^ , pP 

^Fs{eP''') -^1 - 



BERNSTEIN POLYNOMIALS, BERGMAN KERNELS AND TORIC K ABLER VARIETIES 19 

Regarding S, we recall that it is the scalar curvature of the metric gii associated to the 
Kahler form LUps = ^dd{l + l-zp), thus 

S = log(l + \z\Y^ = 2Trgn = 2. 
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